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Abstract 

We make an attempt at unveiling the thermodynamic "signature" of a specific class of 
electronic systems, the so called nearly-flat-band paramagnets and ferromagnets that can 
theoretically be described by appropriate versions of the Hubbard model. 

1 Introduction 

A nearly-flat-band system is an electronic system whose one or more bands in the single-particle 
spectrum can be made arbitrarily narrow by continuously adjusting some parameters of the system, 
without destroying its structure. To make this phrase more explicit, let us consider a tight-binding 
description of an electronic system. In this description a single-particle system is given by a graph 
that consists of a set of sites, with assigned values of external potentials, and a set of bonds 
connecting the sites, associated to non- vanishing hopping intensities. Such a graph is connected if 
any two sites are connected by a sequence of bonds. 

Graphs of nearly-flat-band systems constitute a class of connected graphs, for which it is 
necessary to tune the hopping intensities and on-site external potentials, or only one of those sets 
of parameters, to get one or more degenerate energy levels, the so called non- dispersive bands 
or flat bands, whose degeneracy is proportional to the number of lattice sites. As a result of 
the tuning, the hopping intensities and on-site external potentials satisfy some relations, which 
constitute sufficient conditions for fiat bands. Note, however, that there are connected graphs 
(bipartite graphs) whose spectrum does contain a fiat band, and no tuning of hopping intensities 
is required [1]. 

To each site of a graph there corresponds its neighborhood, that is the set of all the sites 
connected with it by bonds. The topology of a graph is specified by the set of all the neighborhoods. 
It is the topology, and not the geometry, of the underlying connected graph, that decides whether 
flat bands can appear in a system. To the best of our knowledge, there is no characterization 
of the topology of graphs that admit flat bands. Only examples, sometimes of classes, of graphs 
where nearly flat bands do appear are known [21 [3l HI [5l E] . 

By perturbing continuously the sufficient conditions for the flat band, that have to be satisfied 
by hopping intensities and on-site external potentials, we transform a flat-band system into a 
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nearly-flat-band system. This definition of a nearly-flat-band system, via a flat-band system 
is a theoretical, mathematical one. The theoretical limit of non- dispersive band constitutes a 
nonphysical system that cannot be realized in experiment. The residual entropy of multi-electron 
free flat-band system, whose flat band is partially filled, is finite, which violates the Ilird law 
of thermodynamics, and there is no Fermi surface. The Ilird law is violated also in interacting 
systems; see [5l [6] for examples, where the ground-state degeneracy in an interacting case was 
calculated exactly. Moreover, the density of states of those systems is singular. In experiments, 
we can only attempt at constructing a nearly-flat-band system, where by suitably adjusting the 
parameters of the underlying system we can make one or more bands arbitrarily narrow, at least 
in principle. 

Nevertheless, flat-band systems are worth of theoretical studies, since they are more simple for 
analysis than nearly-flat-band systems, and moreover it might happen that some of their features 
are stable against perturbations that transform a flat band into a nearly-flat band. A prominent 
example of such a feature is ferromagnetism, and quite naturally the first model studied was the 
paradigmatic Hubbard model: 

H = 4,aCj,a+J^ni^ani,_„, (1) 

where the sums are over all the sites i,j of the underlying graph, and over projections, a, of 
the electron spin on some axis; tij - the matrix elements of a single-particle Hamiltonian between 
states localized at sites i and j give the hopping intensities and on-site external potentials; c^^, Ci^a, 
stand for the electron creation and annihilation operators, respectively; the term proportional to 
U > represents a strongly screened Coulomb repulsion. When the underlying graph is bipartite, 
Lieb [1] proved the existence of unsaturated ferromagnetism in the ground state (that is, for given 
number of electrons the total spin of the ground state is a fraction of the maximal one), when it is 
a line graph - a proof of saturated ferromagnetism in the ground state (the total spin of the ground 
state is maximal) was given by Mielke [2] , and when it belongs to Tasaki class - the corresponding 
result was obtained by Tasaki [31 H] . In all those cases there is a flat band in the single-electron 
spectrum, and on switching on the Hubbard repulsion the paramagnetic ground state of a free 
multi-electron system turns into a ferromagnetic one, for a special value of electron density (or 
a narrow interval of densities) and any nonzero value of the Hubbard on-site repulsion U. The 
latter statement means that paramagnetic ground state turns into a ferromagnetic one without 
any competition between the kinetic and potential energies (the Hubbard repulsion is needed only 
to lift the macroscopic degeneracy), which is another non-physical feature of flat-band systems. 
The flat-band ferromagnetism discovered by Mielke and Tasaki appeared to be robust against 
perturbations of the flat band. The proofs of this fact for nearly-flat band systems can be found 
in [3 [8]. 

While there is a number of theoretical examples of nearly-flat-band systems (see the papers 
quoted above), we do not know of any measurements performed on real nearly- flat-band systems. 
There have been a few proposals of experimental realizations of nearly-flat-band systems such as: 
atomic quantum wires [H], quantum-dot super-lattices [TU] or organic polymers [H]. However, 
the most promising seems to be a realization as cold atoms in optical lattices [12]. Due to a 
very good control of system parameters in the latter case, such a realization would open new 
possibilities of investigating the mechanism of nearly-flat-band ferromagnetism, inaccessible in 
other experimental realizations, and beyond the scope of present-day theoretical methods. In 
this perspective, it is already interesting and challenging to determine theoretically characteristic 
low-temperature thermodynamic properties of the aforementioned systems. 
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2 The models and the goal 



In this paper, we provide a resume of our investigations of the question: how to recognize a nearly- 
flat-band ferromagnet by means of thermodynamic measurements? In view of the high sensitivity 
of the considered systems to details of the underlying graph, it would be naive to expect that some 
general answer, good for any type of such systems, can be given. We have chosen to concentrate on 
nearly-fiat-band systems described by the Hubbard Hamiltonian, with the lowest band nearly-flat 
and separated from the upper bands by a gap, which we call the principal gap, to distinguish it 
from any other gap in the spectrum of the system. Moreover, we set the Hubbard U to be small 
compared to the value of the principal gap, and the electron density not to exceed ri/b, which is 
the density of those electrons that can be accommodated in the nearly-flat band. For instance, 
the Tasaki models [7] and some related models [13] belong to this class. 

One can argue, however, that our thermodynamic results, to be presented below, hold as well 
for a wider class of systems. This class includes systems where there are two groups of bands: 
low-energy bands and high-energy bands with the nearly-flat band being the highest band in 
the group of the low-energy bands, and separated by a principal gap from the upper bands [9] . 
Another group of systems that belongs to this class are the systems where the flat band sticks 
to a dispersive one in a few points of the Brillouin zone, as in the Mielke class of models (see a 
remark in [T3] and numerical results in [6]). 

To minimize the burden of large volume computer work we have performed calculations for the 
Hubbard model whose graph is a one-dimensional lattice, decorated with additional sites located 
in the middle between the sites of the lattice (known also as a A-chain [13] or a sawtooth chain 
[5]). For appropriate hopping intensities and external potentials the lower band of the model is flat 
and the ground state is ferromagnetic for any ?7 > 0, provided the electron density is nf^jl. This 
model belongs to Tasaki class of models [4] . There are many ways of perturbing the system to get 
a nearly-flat-band system, whose ground state is ferromagnetic only for sufficiently large V . In 
our calculations we chose the Tasaki perturbation [7], whose advantage is that the ferromagnetic 
ground state and its energy are known explicitly. The real matrix elements tjj of the perturbed 
model can be chosen as follows (with the lattice constant of the one-dimensional lattice set to 
unity): = 1, ti,i+i/2 = 1 + s, ti-i/2,i+i/2 = -s, U^i = 2 - s, ti+i/2,i+i/2 = 1 - 2s, for an integer 
i, where s > is the parameter of Tasaki perturbation; up to a factor it amounts to the width of 
the nearly-flat band. One flnds that the width of the lower band is (5_ = 4s, that of the upper 
band - 5+ = 4, and the gap between the bands is e = 1 + s. 

For comparison, we consider also a flctitious noninteracting electron system, not born by a 
Hamiltonian, whose single-particle spectrum consists of two bands, with some dispersion relations, 
separated by a gap. The advantage of this model is that the width of the lower band, 5_, the 
upper band 5+, and the gap between the bands, e, can be varied independently, what facilitates 
observing their impact on thermodynamic properties of the model. 

3 The isochoric heat capacity 

We start our considerations of thermodynamic quantities with the entropy per particle, s(T, v\ as 
a function of temperature, T (throughout the paper T is measured in energy units), and volume 
per particle, v - the inverse of electron density, n. Practically, since this entropy is not directly 
observable, we consider a simply related and accessible to direct measurements quantity - the 
isochoric heat capacity per particle (briefly speciflc heat), cy(T), for specifled values of electron 
density n. The heat capacity at constant volume, is ideally suited for our purposes, since it 
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describes a response of a system to heating, exclusively due to thermal excitations. In contrast 
to other heat capacities, the isochoric heat capacity contains no contribution of a mechanical or 
chemical work done by the system. Such an extra contributions blur the response of the system 
to heating, making it more difficult to identify the system as a nearly-flat-band system. 

Taking into account the variety of graphs and spaces of parameters of nearly-flat-band systems, 
we can claim that ground states of a nearly-flat-band electronic systems are typically paramagnetic, 
that is the total spin Stot = o{Ne), where N^, is the number of electrons. This is true even if we 
restrict the variety of graphs to those that admit ferromagnetic ground states, that is states whose 
Stot = 0{Ne), for suitable values of system parameters. The reason is that if the underlying graph 
admits ferromagnetism, rather limited values of electron densities, nearly-flat band widths and 
Hubbard repulsion U are required to make the ground state ferromagnetic. Quite generally, there 
is no ferromagnetism, in the above sense, if the electron density and/or Hubbard repulsion U are 
too small, and/or the nearly-flat band - too wide. 

A principal gap in a single-particle spectrum induces a gap in the spectrum of the corre- 
sponding many-electron noninteracting system and interacting system described by the Hubbard 
Hamiltonian. Consequently, the spectra of these many-electron systems can be split into a lower 
part and an upper part. As a result, plots of cv{T) consist of two, low-temperature and high- 
temperature, humps. As long as our system is paramagnetic, the low-temperature hump is due 
to low-energy excitations, that is excitations whose energies do not exceed the width of the lower 
part of energy levels (the width of the nearly-flat band in the noninteracting case). The high- 
temperature hump is, in turn, due to high-energy excitations, whose energy does exceed the gap 
between the low-energy group and the high-energy group of energy levels (the principal- gap width 
in a free nearly-flat-band system). Consequently, the position of the low-temperature hump is of 
the order of the width of the nearly-flat band, while that of the high-temperature one is shifted 
by an energy of the order of the principal gap, at least for not too large U. The degree of overlap 
of the two humps, i.e the extent and the depth of the well between them, depends mainly on the 
ratio of the principal gap and the width of the lower band. The larger this ratio is the better the 
separation of the humps. For sufficiently large gaps, the bottom of this well reaches zero and is 
flat. Of course, cv(T) tends to zero if T — )■ or T — )• oo. These features are well illustrated in 
Fig. [H where results for the two-band noninteracting system are shown. In particular it is clear 
that this morphology of the cy(T) plot is no characteristic of a nearly-flat-band system. It can be 
observed in any system of the class described above, interacting or noninteracting, where instead 
of a nearly-flat band there is just a narrow band. 

In contrast to a narrow-band system, in a nearly-flat-band one we can exploit the possibility of 
shrinking the lower band, while preserving the structure of our system. On decreasing 6-, starting 
from sufficiently small value, the low-temperature hump moves, of course, towards zero and almost 
linearly in 6^, its half-width shrinks, while its maximum remains essentially unchanged. These 
effects can be seen in Fig. [T] for the noninteracting two-band model and in Fig. |2]for the particular 
Hubbard model - Tasaki model, defined in previous section. In the theoretical limit of flat band, 
we are left with one, high-temperature hump, which may be separated from zero temperature by 
a visible plateau of zero value, provided the principal gap is large enough. This is because there 
is practically no low-temperature excitations from the flat band. 

Now, suppose that the underlying graph, the electron density, and the Hubbard repulsion U 
are such that for those widths of the nearly-flat band that are smaller than some threshold value, 
6p-f (which depends on U), the paramagnetic ground state changes into a ferromagnetic one. 
Starting measurements with the values of (5_ somewhat larger than 5p-/, and then repeating them 
for a decreasing sequence of values above one observes the above described "evolution" of the 
low-temperature hump (see Fig. [2]). However, below 6p^f the gapless low-temperature excitations 
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Figure 1: Specific heat versus temperature, cv(T), of tlie noninteracting two-band model, witli 
5+ = 1 and 5 = 6 (in arbitrary energy units), and partially filled lower band, for a decreasing 
sequence of lower-band widths 6- and for the flat band. Temperature is in a logarithmic scale. 

of a paramagnet (of hole-particle type) turn into ferromagnetic excitations - the magnons. This 
transition results in pinning the low-temperature hump at the temperature of the order of the 
width of the magnons band, which is determined by U. Further decrease of down to (5_ = 
0, i.e. the limit of flat-band ferromagnet, does not bring any significant changes to the low- 
temperature hump, neither to its position nor to its shape. In Fig. [3] we see a transient region 
between the paramagnet and the ferromagnet, with the values of 5_ greater than 5p-/, but smaller 
than the values of 6- for which a typical for a paramagnet, linear in 6-, "motion" of the low- 
temperature hump towards zero is observed. In contradistinction to the overall shape of cv{T) 
plots in paramagnetic and ferromagnetic states, we expect that the size (measured by a range of 
(5_) of the transient region and the shape of cy(T) plots in this region are rather sensitive to the size 
of the system used for calculations. Finally, in Fig. HJ we see the pinning of the low-temperature 
hump for 5_ smaller than 5p-/. 

4 Heat capacity at constant chemical potential 

In the previous section we discussed the heat capacity per particle at constant volume, which is 
proportional to the second derivative with respect to temperature, T, of the fundamental ther- 
modynamic function, the Helmholtz free energy per particle as a function of T and volume per 
particle, v. An analog of this thermodynamic quantity, biased by finite-size effects, can be calcu- 
lated numerically for small Hubbard systems, with fixed volume and particle number, by means of 
the canonical ensemble. Formally, the heat capacity at constant chemical potential, per unit vol- 
ume, c^, is another quantity of this kind. It is proportional to the second derivative with respect to 
temperature of the fundamental thermodynamic function, the grand-canonical potential per unit 
volume (i.e. minus the pressure) as a function of temperature and chemical potential, /i. This is, 
however, an unusual heat capacity, since it refers to open systems that are in thermal equilibrium 
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Figure 2: Specific heat versus temperature, cy(T), of tlie particular Hubbard model - Tasaki 
model, defined in the paper, that consists of ten sites with periodic boundary conditions, and five 
electrons (quarter filling), for t/ = 0.1 and for the range of lower-band widths 5_, where a typical 
paramagnetic behavior is observed. Temperature is in a logarithmic scale. 



with particle reservoir of given chemical potential. Consequently, in a plot of versus T, different 
points correspond, in general (depending on /i, the electron density at constant /i, n^(T), can be 
monotonic or not), to systems with different amount of matter. An analog of this thermodynamic 
quantity, biased by finite-size effects, can be calculated numerically for small Hubbard systems, 
with fixed volume and chemical potential, by means of the grand-canonical ensemble. Of course, 
c^(T) can be related to ncy(T), which is the isochoric heat capacity per unit volume, by the 
thermodynamic identity: 

CAT) = „,(r)c.(T..(r.,)) + r (^)' (^)" ■ P) 

The difference between those two heat capacities, — n^cy, is born by the chemical work that 
an open system does exchanging matter with matter reservoir when heated. 

In a free system, whose spectrum consists of a fiat band only, cy(T) vanishes identically, since 
there are no excitations in this system. In contradistinction to cy(T), c^{T) can be nonzero, due 
to the chemical work, the heated system does exchanging matter with a matter reservoir. In such 
a system the electron density is a function of the activity ^ = exp(/i/T) only; we used here the 
fact that we can always set the energy of the fiat band to zero. Hence, the second term of identity 
([2]) - the chemical work term can be written as 

ein^(O^>0. (3) 

In a free system, with only a flat band, dn{S,)/d^ > 0, therefore the chemical work term, is nonzero 
unless the chemical potential coincides with the energy of the flat band. Consequently, there is 
just one hump in a plot of c^{T) for /i 7^ 0, located at low temperatures, if \fj,\ is small. The size 
and the position of this hump is sensitive to the value of the chemical potential. 
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Figure 3: Specific heat versus temperature, cy(T), of tlie particular Hubbard model - Tasaki 
model, defined in the paper, that consists of ten sites with periodic boundary conditions, and five 
electrons (quarter filling), for U = 0.1 and for the range of lower-band widths 5_, greater than 
6p-f but smaller than the values of 6- in Fig. [2l Temperature is in a logarithmic scale. 

If a fiat band is accompanied by some upper bands, separated by a gap from the fiat one 
(like in the systems considered in the previous section), the above statements concerning low- 
temperature hump remain qualitatively true at sufficiently low temperatures and for sufficiently 
small Additionally, such a system contributes significantly to c^{T) at high-temperatures; 
there is a high-temperature hump in the plot of c^{T), which coincides essentially with ncv{T) 
- the isochoric heat capacity per unit volume. Small guarantees that the low- and high- 
temperature humps are well separated. 

A weak perturbation that turns a flat-band system into a nearly-fiat-band one, could have 
influenced significantly only the low-temperature hump. Unlike the low-temperature hump in 
cv{T) plot of a free nearly-fiat-band system, which is due to thermal excitations from a nearly- 
flat-band, and therefore is sensitive to the width of this band (its position moves towards zero 
temperature as the perturbation decreases), the low-temperature hump in c^(T) is due to exchange 
of matter, and its position and shape are quite insensitive to the width of the nearly-flat band. 
Thus, a weak perturbation of the flat band does not change essentially the two-hump c^(T) plot 
of the flat-band case. The above observations are well illustrated in Figs. O |6l [3 

Now, consider the Hubbard model ([T]) whose graph is specified in Section 2. In the absence of 
perturbation (s = 0), the lower band in the single-particle spectrum is flat and is separated by a 
gap from the upper band. There is a basis of the flat-band eigensubspace that consists of localized 
eigenstates. In terms of Slater determinants of those localized single-particle eigenstates, one 
can construct a basis of the ground-state eigensubspace of the Hubbard model, for any Hubbard 
repulsion U > 0, provided the number of electrons does not exceed half-filling of the flat-band 
Derzhko et al [S] have demonstrated that those bases, as well as ground-state-eigensubspace 
bases of some similar flat-band Hubbard models [6], can be mapped onto a fictitious hard-core 
lattice gas. The ground-state is paramagnetic, and the many-electron excited states are separated 
by a gap from the ground state, for any filling that is smaller than half-filling of the flat-band. 
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Figure 4: Specific heat versus temperature, cy(T), of tlie particular Hubbard model - Tasaki 
model, defined in the paper, that consists of ten sites with periodic boundary conditions, and five 
electrons (quarter filling), for t/ = 0.1 and for the range of lower-band widths 5- that are smaller 
than and for the flat band. Temperature is in a logarithmic scale. 

Consequently, for such fillings, at sufficiently low temperatures that make excitations above the 
ground state very unprobable, the response of the system to heating amounts essentially to that 
of the hard-core gas corresponding to the macroscopically-degenerate ground state. Since there 
are no excitations in a closed hard-core gas, cy vanishes identically. In contrast, the chemical 
work term dH]), which in the case under consideration and in many other cases can be calculated 
analytically, is nonzero for 7^ 0. Therefore, there is just one hump in a plot of c^{T) for /i 7^ 0, 
at low temperatures. 

A weak perturbation of the fiat-band (s ^ 1) does not bring essential changes to this low- 
temperature hump. It is given essentially by the chemical work term of the hard-core gas corre- 
sponding to the ground state of the fiat-band limit of the considered nearly-fiat-band Hubbard 
model. Therefore, all the above described features of the low-temperature hump of c^(T) plot in 
free nearly-flat-band systems, in particular its insensitivity to the width of the nearly-flat band, 
remain valid in the considered Hubbard model. These observations apply also to other Hubbard 
systems mentioned in Section 2. It is the derivative dn{^)/dC, that differentiates between hard-core 
gases and/or flat band systems and makes the shape of the low-temperature hump specific for a 
system. However, those differences are not dramatic; they are hidden in fine details of the shape 
of the low-temperature hump (see the plots in [B]). 

Naturally, at high temperatures, there develops a high-temperature hump due to excitations 
above the ground state, much like in the case of a free systems considered above, and it amounts 
essentially to the specific heat per unit volume, ncy(T). The overall picture of c^(T) plot in the 
considered Hubbard models is qualitatively the same as those in free nearly-flat-band systems con- 
sidered in the previous paragraph. The described above features of c^{T) plot are well illustrated 
by the plots displayed in [5], for the Hubbard model considered in this paper, and in [6] for other 
Hubbard models with flat or nearly-flat-bands. The data for all those plots were obtained from 
exact diagonalization of small Hubbard systems, with the number of sites in their graphs between 
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Figure 5: Plot of c^{T) for the noninteracting two-band model, with 5+ = 1 and £ = 6 (in arbitrary 
energy units), with the flat-band energy and the upper edge of the corresponding nearly- flat band 
set to zero. Temperature is in a logarithmic scale. 

10 and 20. Apparently, looking at a c^(T) plot on can hardly infer whether the underlying system 
is a nearly-fiat-band (fiat-band) Hubbard system. 

5 Summary 

Summing up, we propose an answer to the question asked in the title of our report. We argue 
that the isochoric heat capacity per particle, cy(T), is a good candidate, a sufficiently sensitive 
thermmodynamic quantity, to measure. At the heart of our answer is the fact that in the case 
of nearly-fiat-band systems we can vary the width of the nearly-flat lower band, and make it as 
narrow as we wish. Then, we can watch how cy(T) plot "evolves" with the decreasing width of the 
lower-band, at sufficiently low temperatures. This low-temperature "evolution" provides a signa- 
ture of a nearly- fiat band paramagnet, paramagnet-ferromagnet transition, and nearly-fiat-band 
ferromagnet. A thermodynamic measurement performed only for a single value of a nearly- flat- 
band width is not sufficient for this purpose. We demonstrate also that the heat capacity at 
constant chemical potential, per unit volume, versus T, c^(T), is not a suitable quantity to mea- 
sure, since it depends weakly on the width of a nearly-flat-band and is sensitive to the value of the 
chemical potential kept constant. The shape of c^(r) plot is no characteristic of a nearly- flat-band 
(flat-band) Hubbard system. One can ask naturally, whether such thermodynamic characteristics 
of closed systems like the coefficient of thermal expansion or the isothermal compressibility can 
be used to recognize nearly-flat-band paramagnets and ferromagnets. Our studies of the nearly- 
flat-band two-band model suggest that those quantities can be used to detect nearly-flat-band 
Hubbard paramagnets. However, those quantities are not well defined and cannot be calculated 
for small lattice systems. Therefore, we have no data for Hubbard systems with ferromagnetic 
ground state to check if those quantities are suitable also for detecting nearly-flat-band ferromag- 
nets. 
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Figure 6: Plot of nfj,(T)cv(T,v(T, ^)) versus T, for the noninteracting two-band model, with 
= 1 and e = 6 (in arbitrary energy units), with the flat-band energy and the upper edge of the 
corresponding nearly-flat band set to zero. The dashed line coincides with the continuous one. 
Temperature is in a logarithmic scale. 
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